Abstract 



We analyze the thermodynamical behavior of black holes in closed 
finite boxes. First the black hole mass evolution is analyzed in an 
initially empty box. Using the conservation of the energy and the 
Hawking evaporation flux, we deduce a minimal volume above which 
one black hole can loss all of its mass to the box, a result which agrees 
with the previous analysis made by Page. We then obtain analo- 
gous results using a box initially containing radiation, allowed to be 
absorbed by the black hole. The equilibrium times and masses are 
evaluated and their behavior discussed to highlight some interesting 
features arising. These results are generalized to black holes + 
thermal radiation. Using physically simple arguments, we prove that 
these black holes achieve the same equilibrium masses (even that the 
initial masses were different). The entropy of the system is used to 
obtain the dependence of the equilibrium mass on the box volume, 
number of black holes and the initial radiation. The equilibrium mass 
is shown to be proportional to a positive power law of the effective vol- 
ume (contrary to naive expectations), a result explained in terms of 
the detailed features of the system. The effect of the reflection of the 
radiation on the box walls which comes back into the black hole is ex- 
plicitly considered. All these results (some of them counter-intuitive) 
may be useful to formulate alternative problems in thermodynamic 
courses for graduate and advanced undergraduate students. A hand- 
ful of them are suggested in the Appendix. 
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In a seminal work of the 70s, Hawking [1] , showed that black holes are capa- 
ble of emitting radiation, and evaluated the quantum process by which these 
black holes lose mass due to vacuum polarization (induced by the gravita- 
tional energy of the black hole). This important work and follow-up contri- 
butions, together with Bekenstein's arguments [2] about the entropy of black 
holes, prompted a new discipline of black hole research: thermodynamics of 
black holes. Today, the study of black hole thermodynamics is very active 
and several arguments indicate that it must embrace the partial (or perhaps 
complete) marriage of quantum mechanics, thermodynamics and gravitation. 

The original Hawking analysis showed that the particles created display 
a thermal spectrum f{E,T) = -^pw^ with (3 — l/ksT, and a temperature 
determined by the black hole mass M given by 



see [3] for details. Since this spectrum is thermal, the luminosity of the 
black hole can be calculated using the Stefan-Boltzmann law, using L oc 

T^r-g^, where the gravitational radius r-g oc M squared plays the role of the 
emitting area. Using these definitions, L{M) oc is readily obtained. 
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In addition to "purely academic" cases (e.g. black holes inside ideal 
boxes), these emission (absorption) relations are in principle useful for inves- 
tigating the properties of astrophysical and/or primordial black holes (hear- 
after PBHs) and their cosmological consequences. For example, some short- 
timescale gamma-ray burst events have been possibly explained by evapo- 
rating PBHs in a terminal (explosive) phase. In fact, the derived relations 
between period, mass and spectra are consistent with the PBH model, see 
[4]. 

Even in the simplest cases (the black holes interacting with ambient radia- 
tion) it is interesting to understand what determines the equilibrium features 
of the system and how is it achieved. We shall perform an analysis restricted 
to finite volumes with the aim of illustrating some novel features of the ther- 
modynamics of black holes using gedanken experiments. These excercises 
serve also as a prchidc to the more complicated situation, in which evapora- 
tion of PBHs must be studied in an expanding universe. However, and quite 
independently of these further advanced considerations, "boxed" black holes 
display interesting features which need careful clarification. 



2 Classical absorption in an initially cold box 

Let us consider the following gedanken experiment: we start with a black 
hole and put it within a finite box (with linear size L) . In the beginning, the 
box has no radiation, i.e. its temperature is zero. Then, since the black hole 
is hotter than the environment, it will emit thermal radiation, and thus this 
closed box will be increasingly filled with this radiation shortly afterwards. 

The evolution of the mass of the black hole is described by the well-known 
expression which balances the energy loses to the Hawking luminosity oc 

dM^_C_ 

dt ^ ^ 

with an initial condition of the box Tj.ad{t — U) — 0. The constant C 
depends on the degrees of freedom allowed to be radiated by the black hole 
[1] and is set to C ~ 10^^(?^s^^ throughout this work, and the initial mass 
is Mj. We assume that the box linear size is much larger than the black 
hole gravitational radius L ^ to avoid a complicated non-linear feedback 
between the geometry and the radiation field, which would lead us far beyond 
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the purely thermodynamical approach. 

Since the box is closed, the total energy is conserved. Therefore it is 
easy determine a relation for the black hole mass and the thermal radiation 
contained within the box from 



Qrad{Trad)Vbo. + M{t) = M, (3) 

where t > ti (we use natural units to set the speed of light c equal to 1). 
The radiation density is given by QradiT) = a^T'^ with a* ~ 8 x lQ^^^gcm~^K~'^. 
After some time, the box temperature is 



Trad{t) 



Mi - M{t) 



nV4 



(4) 



At an even later time t > ti, eq.(2) should not simply contain the Hawk- 
ing term, but an absorption term is also needed. The simplest form of term 
is generally constructed as the product of the incoming flux times the gravi- 
tational cross-section of the radiation faUing into the hole, Bgrad{T)M'^, with 
B — ^^^^-{or^\ in natural units, see Ref. [5] for details and discussion). 
Therefore the complete differential equation for the mass is given by 

'f^-^ + B,UT)M^ (5) 

From this equation, it is clear that thermodynamical equilibrium between 
the black hole plus and the ambient radiation will be given by M = 0, and 
thus the mass of the black hole in equilibrium is 

= M,{te,) = ^ (6) 

radKj'eq) 

where D — [C/a^B)^^'^ ~ 2 x IQ'^^gK. Solving the set of equations above 
we track the black hole mass evolution as it approaches the equilibrium (and 
determine the conditions for this equilibrium). 

If we substitute the eq. (4) into eq. (5) we obtain 

The evolution of the mass of eq.(7) will be given by solving the integral 
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J Mi -Ki + + KsM^ ~ ^ ' 

and the constants are Ki — C, K2 — BMi/Vbox and — —B/V^ox- 

In the hmit V^ox » fg^ the solution of eq.(7) approaches M{t) ~ Mj[l — {t/tevap)t^^ 

with tevap — 'it'' ^ expected. 

The numerical solution of the integral cq.(7) is shown in Fig.l for different 
box sizes but the same initial black hole mass Mj. In order to obtain the 
conditions for the equilibrium we solve the algebraic equation Trad{teq) = 
THawiMeq), rewritten as 



Mi - M{t 



eqj 



a*Vb 



box 



1/4 ^ 



M{t 



(9) 



eqJ 



Inserting the respective constants above, eq.(9) becomes 



[1 - /i(te,)]/i'(te,) = D{Vhox/cm^) X (Mq/M,)' (10) 

where /x(teg) = and D - 4.2 x IQ-^^. 

The interpretation of the solution is quite evident, when the volume of 
the box is such that the right side of eq.(9) is larger than one, the black hole 
evaporates completely in the absence of additional quantum corrections at 
the Planck scale. But if D{Vbox/cm^) x (Mq/Mj)^ ~ 1, the size of the box is 
small enough to stop the evaporation, and the system black hole+radiation 
achieves thermodynamical equilibrium before its mass vanishes completely. 

The hnear size of the box below which the black hole achieves its equihb- 
rium is given by 

Lc(Mi) - 1.3 X lO^^iMi/M^f^^ cm (11) 

The time for reaching equilibrium teg is determined by the initial mass 
and the size of the box only, i.e. t^g = teq{Mi, Vbox)- For the same initial mass, 
smaller boxes will contain black holes with larger final masses at equilibrium 
(see Fig. 1). Our analysis agrees with the previous work by Page [7]. From 
now on, we define the critical volume by V^it(M) = Lc{M)^ using eq.(lO), 
with its clear physical interpretation given above. In the next sections we 
revisit the derivation of eq.(lO), extend it to more complicated cases and 
consider causality features. 
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Figure 1: The approach to equihbrium of black holes in finite boxes. The 
curves represent qualitatively the temporal behavior of the mass of the black 
hole for two different values of the box size V^oa; (initially devoid of radiation) ; 
with the initial value of the mass Mj held fixed. 

3 Classical absorption in a closed box with 
initial radiation 

Let us now evaluate the behavior of the black hole mass when introduced in 
a box with a non-zero initial radiation content. We expect the black hole to 
achieve thermodynamical equilibrium earlier than in the case without initial 
radiation, and with a smaller equilibrium mass (considering boxes with the 
same size). In this case, is easy show that 

J Mi -Ki + KiM^{Ei - M) ~^ ^ ' 

where E'j = Mj + a^Ti'^Vbox and = B/Vbox- The plot of M{t) as a 
function of time is actually similar to the former case. 

As before, we impose the equilibrium condition in to obtain the relation 
between the black hole mass and the box volume. Following the same steps 
as before, we obtain 

/x*(te,)'[F(T„ M,) - /i*(te,)] = D{Vho./cm^) X (Mo/M,)' (13) 

where F(Ti,Mi) — and the teq is different from the previous case. 
Actually the black hole achieves thermodynamical equilibrium sooner, as 
expected. 

4 N black holes plus radiation in a closed 
box 

The generalization to black holes seems quite straightforward, although it 
will become clear that the initial states and other details must be carefully 
defined to achieve consistent results. First, we shall consider just two black 
holes immersed in a closed box (with constant volume) in the following initial 
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situation: one black hole has initial mass Mi and the other black hole is more 
massive than the first, M2 > Mi. Initially the box has no radiation, and 
therefore these objects begin to evaporate immediately. Moreover, we shall 
consider that V^ox ~ Kr-it(-^i)- Then at the initial time ti we have 



dt Ju Ml 

dMi\ C 



dt Ju Ml 



(15) 



Since Mi < M2 one may consider |(Mi)emp| ^ |(-^2)emp| initially. Af- 
terwards, when the box starts to be filled with the emitted radiation, some 
energy can be absorbed by the black holes. Therefore, their evolution will 
be given by 

^ +BgUT)M,' (16) 



dt M| 

dMi _ C 
IT ~ ~Mf 



+ Bgrad{T)Mi^ (17) 



Mi{t) + M2{t) + VboxQradiT) = Mi{ti) + M2{ti) (18) 

where the last condition displays the conserved energy of the box. Using 
this constraint, we can describe this system by the following set of equations 

combining eqs.(17) and (18) and using the conservation of energy yields 



-C 



1 1 

M? ^ M| 



+ BQradiMl + M|] + VboxQrad = (21) 



where = Mi + M2 is the initial energy. 

Note that the set of coupled equations does not include the effect of 
black hole motion due to their mutual gravity. Motion is hkely to affect the 
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emission/absorption properties of the black holes in the relativistic regime 
[8], [9]. Thus, the analysis is strictly valid whenever Vbh <^ c. 

Let us take a look at the global behavior of the solutions. First, it is easy 
to show that the equilibrium between two black holes is possible. If we impose 
that the thermodynamical equilibrium will be achieved, then asymptotically 

dm^dj^^ (22) 

From these equations ^Vbox = — Mi — M2] immediately follows, 

and the same is true for the other black hole. Therefore, it follows that 

Mi,eg = M2,eg (23) 

The same results of these equilibrium masses could have been obtained 
from the expression of the entropy. In fact, more complicated cases can be 
worked out either using the above approach or using that the entropy must 
be an extreme for a system in equilibrium. For instance, let us consider N 
black holes enclosed in the box. The total entropy for this system is given 

by 

Stotaiit) = \j2 Mt) + g*TUt)Vi,o. (24) 

i 

Since the horizon area for the i-th black hole is Ai — 'inVg/ we have 

An 

Stotaiit) = TT4 E Miitf + g.T'radimox (25) 

In equihbrium T^ad = ^b/ii = Ti,^2 = ... = Tfe^i = 8^;rp^> and therefore 

S^ax{N, Vbox) = /9*iV/ieQ^ + h/l^eq^{Vbox/cm^) (26) 

where f3^ = 2.5 x 10^°, = 8 x lO^^a* and i^eq = {Meq/lO^^ g) the 
equilibrium mass scaled to a convenient reference value. 

Extremizing the entropy ^^§(f^^ = and solving for Meg yields 

Meq{N, Hox) ~ 5.4 X 10'^a,'/'N-'/\V,oJcm^f' g (27) 
with a* = 26.7g^. After substituting back into S^ax we find 
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Figure 2: Graphical solution of the non-linear equation for the equilibrium 
mass (eq.30) as a function of time for a fixed value of Vbox and a different 
number of black holes A^. It may be said that the equilibrium mass "feels" 
the presence of other black holes. 



SmaAN: Hox) ~ 1-2 X N^/' {VboJ cm^Y^' (28) 

Since = nVbox by definition, we have Smax ~ 1.2 X lOV^VboJcm^). 
Therefore, the thermodynamic approach leads us to conclude that the equi- 
librium black hole mass is bigger for bigger boxes. This result is quite counter- 
intuitive, but it has a reasonable explanation given below. To confirm the 
dependence first, we resort to another method. 

Let us rewrite the equation for the mass of the i-th black hole as 

^ - -XfW + - ^^^(*)] (29) 

dt Mi{ty Vbox 

where we prepare the system with N black holes with the same initial 
mass (this particular condition is not over-restrictive and simplifies the alge- 
bra) . Prom Mj = we obtain 

/^e/ - -/^e/ = -^{Vbox/cm^) (30) 

where K ~ 3.5 x 10^^. This algebraic equation has non-trivial solutions 
which can be found numerically. For the sake of definiteness we adopt = 1, 
— 10, and a small Vbox = lO^cm^, which yields fj,eq ~ 0.137. If wc enlarge 
the box, say Vbox = lO^cm^ instead, then ~ 0.244, and so on. If we fix 
the size of the box instead, say Vmx = 2.8 x lO^cm"^, and change the number 
of black holes, the equilibrium mass will refiect the effect of the presence 
of other black holes accordingly. For example, taking Hi — 10 and N — 1 
wc have /deq ~ 1-027 but for N = 10^ we obtain //e, ~ 0.178 (Fig. 2). 
These results agree with the previous entropy analysis, and confirm that the 
coupled equations correctly describe the evolution of the system. It is quite 
interesting to propose the construction of a plot showing the behavior of the 
solutions and a discussion about the reasons for that behavior. 
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As a further related work it is proposed to show that if we consider the 
existence of an initial thermal radiation with energy Ei, the effective volume 
becomes 

VefAN, Mi) = (31) 

which holds for N black holes with the same initial masses. This effective 
volume must be used instead of (Hox/-^) in eq.(26) and Meq{N, Vbox) becomes 
^ r v^o. 

The puzzling growth of the equilibrium mass in these systems must be 
further explained in more detail to improve the understanding. To proceed, 
let us consider two boxes with volumes Vhox,2 and Vhox,i where the second 
box is larger than the first. Let us immerse two black holes with the same 
initial masses in both boxes. Therefore, the energy content is forced to be the 
same. Then, in equilibrium we will have (dropping some inessential numerical 
factors) 

Ebox = 2Meg,i + T^ad,l^box,l = 2Meg^2 + T^ad,2^box,2 (32) 

Since the equilibrium mass is given by Mgg oc [Ho-r]^^^ for both boxes, 
then we have 

V'box,!'^^ - ybox,2^^^ = T^ad,2^box,2 " T^ad,l^box,l (33) 

both terms are negative due to the assumption Vbox,2 > Vbox,i then we 
obtain from the r.h.s. of this equation 

^)'<(^)<1 (34) 

then Trad,2 < T^ad,!- Therefore, the bigger box will have smaller radiation 
temperature. Since in equilibrium Trad2 = Tbh2 oc -rr — then the mass of 

' ' -'"eg, 2 

the second black hole must be bigger than the first, in agreement with our 
previous analysis. Previous work by Page [7] (specifically his eq.(lO)) has 
explored this problem, which has been solved here using the conservation of 
energy only, with the absorption terms explicitly present. 
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5 Causality considerations 



The above considerations implicitly assumed that the black hole absorption 

starts instantaneously as the former is created inside the box. One might 
wonder whether a black hole immersed within an initially empty box is actu- 
ally able to absorb ambient radiation immediately; since the Hawking radia- 
tion emitted by it will be available to absorption only after ~ L/cas required 
by causality. 

In order to gain some insight, we use the fact that if the box is large 
enough, when the radiation comes back into the black hole, the object has 
evaporated completely. Then, a critical condition is obtained imposing this 
time interval to be bigger than the timescale for evaporation. Thus 



Then, any box (without initial radiation) containing a black hole with 
mass M whose linear size L is bigger than the critical size Lcrit above would 
let the black hole to evaporate completely, thus precluding equilibrium. Note 
that the functional dependence and numerical value are very different from 
the naive approach used previously where we have considered that the ab- 
sorption term turns on instantaneously (that is, as soon we immersed the 
black hole in the box). If we boldly compare Lcru with the size of the Hubble 
horizon today ~ lO^^cm. we deduce that a PBHs with masses smaller than 
~ 2 X lO^^g can evaporate completely. However, it is clear that in a realistic 
case other sources of energy must be considered and our comment is intended 
just for pedagogical purposes. 

The result of eq.(35) can be interpreted as follows. If we put one hot black 
hole (at ti) with initial mass inside a closed box with volume Vbox- it will 
evolve according to M = —C/M'^ for a time ti < t < ^i^. But from t > 
Lbox/c on, the black hole has to absorb radiation from the box (originated by 
itself!). Then the absorption term oc T^Hox contributes to avoid its complete 
evaporation. If the box is slightly larger (but smaller than the critical volume 
defined by eq.(lO)) the onset of absorption is delayed, then this black hole 
will be hotter than the previous case. This black hole is then filling the box 
with higher temperature radiation, and after some time it will be absorbing 



L 

7 ^ 3C 



(35) 



Therefore, it follows that Lcrit{M) ~ lO^^cm 




11 



radiation at higher temperature, hence the absorption term oc T'^M'^ will be 
larger. This term drives the black hole into the larger equilibrium mass, as 
given by eq.(27). It should be remembered that a maximal number for black 
holes enclosed in a box with volume Vbox must exist, a naive estimate of N^^ax 
can be obtained from the close-packing condition Nj^axf^^^eq) ~ Hox- 

6 Conclusions 

In this work we have analyzed the behavior of black holes+radiation sys- 
tems, considering some gedanken experiments within finite boxes (initially 
with and without radiation) (see [10] for a clear pioneering discussion of 
these issues). These 5'e(ianA;en experiments are important to understand the 
evolution of black holes in more complex situations, since it is hoped that we 
can always approximate the thermodynamical behavior of the universe with 
that of a finite box. In these simple cases the energy available for absorp- 
tion by the black holes is finite at all times and there is no need to consider 
other sources. Strictly speaking, and with an eye on the actual cosmological 
formation and evolution of black holes, it may be argued that the problem 
is oversimplified, since actual primordial black hole masses must be formed 
with a substantial fraction of the horizon mass [6], and therefore finite size 
effects should be considered from the scratch. The actual situation is much 
more complicated but not impossible to tackle. The real difficulty, which can 
be viewed as an advanced exercise dealing with the (generalized) second law 
of thermodynamics not related to laboratory systems, is that (unlike labo- 
ratory gases) the black hole "gas" will have varying masses and therefore a 
complicated kinetic equation must be used to describe their evolution. This 
provides a concrete example, possibly realized in nature, of a variable-mass 
gas needing a deep understanding of statistical mechanics principles for its 
very formulation. The simplest approach taken here was enough to show 
that one black hole plus radiation can achieve thermodynamical equilibrium 
if the box volume is smaller than a critical volume (in agreement with a 
previous treatment made by Page [7]). Some surprises arise in more general 
cases where we have two (or more) black holes (plus radiation) within these 
boxes. 

Even in their simplest versions the cases of black holes in boxes are very 
instructive to analyze and may serve as an good issues for a graduate course 
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of thermodynamics/statistical mechanics (see [11] for a discussion of related 
formal aspects of black hole thermodynamics) as an alternative to tradi- 
tional problems dealing with the approach to equilibrium. They pose several 
questions of deep physical meaning, which are also strongly entangled and 
sometimes puzzling to interpret. 
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7 Appendix 

A complementary set of simple problems intended to reinforce the conceptual 
aspects of black hole evaporation. Except for the problem 4, their mathe- 
matical complexity is quite straightforward. 

1) Using eq.(2), prove that the timescale for evaporation for one black 
hole is proportional to M^. Verify that for one black hole with the solar 
mass (M ~ 2 X 10^^ g); tevap{M) ~ l^f^years. Evaluate this for a initial 
mass ~ 10^^ (? and compare with the age of the universe. 

2) Verify the expression eq.(6) for the critical mass. Evaluate this number 
today, considering that Trad{to) ~ SX. 

3) Derive eq.(ll). Use for the hnear size of the box the present size of 
the universe. Evaluate the corresponding mass at equilibrium. 

4) Deduce the set of eqs.(19)-(21) and interpret the solutions graphically. 

5) Deduce the Mgq{N, Vhox) for the cases without initial radiation and 
with a radiation filling the box. Interpret these results and compare with the 
cosmological situation plugging representative numbers for the latter. 

6) Discuss how all these results may be modified if the box expands 
following a known temporal dependence. 
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